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A snake (of length L) is a (continuous) piecewise (^-curve S : [0, L] — * 
M. d , parameterized by arc-length and whose "tail" is at the origin (5(0) = 0). 
■ Charming a snake consists in having it move in such a way that its "snout" 

/"*\ [ follows a chosen C 1 -curve 7(t). The snake charmer algorithm, initiated in 

|Ha2| for polygonal snakes and developed in |Ro| in the general case, works as 
| follows. The input is a pair (5,7), where: 



(i) S : [0, L] -> R d is a snake of length L, 

(ii) 7 : [0, 1] -> R d is C^curve with 7 (0) = S(L). 



The output will then be a continuous 1-parameter family St of snakes of length 
<^ \ L satisfying So = S and S t {L) = ^(t). This algorithm, described in Section [5] 

CNJ ■ below, is Ehresmannian in nature: the output is a horizontal lifting for some 

connection. A holonomy phenomenon for closed curves 7 then occurs: having 
7(1) = 7(0) does not imply that Si = Sq. Given a snake 5, one can input in 
the snake charmer algorithm the pairs (iS, 7) for all possible C 1 -loops 7 at S(L). 

■ The snakes Si obtained this way form the holonomy orbit of iS. The purpose of 
this note is to study these holonomy orbits, proving that, in good cases, they 
are compact smooth manifolds diffeomorphic to real Stiefel manifolds. 

The paper is organized as follows. After some preliminaries in Section Q 
we give m Section a survey of the snake charmer algorithm and of its main 
^ I properties (more details are to be found in [Rop . Section [3] presents our new 

•rH . results about holonomy orbits with their proofs. Section 0] contains examples. 

/\ ' The first author is grateful to Giorgi Khimshiashvili for a stimulating dis- 

■ cussion during the conference. 



1 Preliminaries 

1.1 Let L be a positive real number and let V = {0 = So < si < • • • < sat_i < 
sat = L} be a finite partition of [0, L\. If ( A, d) is a metric space, a map 
z : [0, L] — > X is said piecewise continuous for V if, for every i = 0, . . . , N — 1, 
the restriction of z to the semi-open interval [si,Si+i) extends to a (unique) 
continuous map Zj defined on the closed interval [sj,Si+i]. In particular, z is 
continuous on the right and the discontinuities, only possible at points of V ', 
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are just jumps. We denote by Cp([0, L],X) the space of maps from [0, L] to X 
which are piecewise continuous for V; this space is endowed with the uniform 
convergence distance 

d(zi,z 2 )= sup {zi(s), z 2 (s)} . 
se[o,L] 

When V is empty, the map z is just continuous andC2([0, L],X) = C°([0, L], X). 
The map z <— > (zi, ■ ■ ■ , Zn) provides a homeomorphism 

iV-l 

C°([0,£],X)« J] C°({ Si , Si+1 },X) . (1) 

i=0 

If X is a Riemannian manifold, the space C°([sj, Si+i], X) naturally inherits 
a Banach manifold structure (see |Ee|1. so Cp([0, L),X) is a Banach manifold 
using JIJ. The tangent space T z Cj,([0, L],X) to Cp([0, L], X) at z is the space of 
those v € Cp([0, L], TX) satisfying pow = z, where p : TX — > X is the natural 
projection. 

1.2 The unit sphere in K d centered at the origin is denoted by S d_1 . Let 
Mob(<i — 1) be the group of Mobius transformations of It is a Lie group 

of dimension d(d + l)/2, with SO(d) as a compact maximal subgroup. Its Lie 
algebra is denoted by m'6b(d — 1). For / d £ M d , we define a 1-parameter 
subgroup 17 of M6b(d - 1) by 

rV — 1 V 

t = fv °Pt o( Pv 

where 

• ip v : S d_1 — * M. d is the stereographic projection sending v/\v\ to 00 and 
-v/\v\ to 0; 

• p v t : R d -> R d is the homothety p v t (x) = e^x. 

Thus, F( is a purely hyperbolic flow with stable fixed point v/\v\ and unstable 
fixed point —v/\v\. We agree that Fj = id. Let C v S mob(<i — 1) such that 
Tj = exp(tC v ) (with Co = 0). The correspondence v i— > C v gives an injective 
linear map \ '■ ~~ * mob(d — 1); its image is a d-dimensional vector subspace 
Ti. of mdb(<i — 1), supplementary to 50(d), and 7i generates m'6b(d — 1) as a Lie 
algebra. Let A w be the right invariant distribution on M6b(d — 1) which is 
equal to TL at the unit element. 

The 1-parameter subgroups may be used to built up a diffeomorphism 
* : R d x SO(d -l)-^> M6b(d - 1) defined by p) = T\ ■ p. 

2 The algorithm 

In this section, we give a survey of the snake charmer algorithm and some of its 
properties. For details, see |Ro| . 
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2.1 Fix a positive real number L and a finite set V C [0,L] as in 11.11 Let 
Conf = Cp([0, L], § d_1 ), with its Banach manifold structure coming from the 
standard Riemannian structure on § d_1 . The inclusion of § d_1 into R d makes 
Conf a submanifold of the Banach space II^lo 1 ^°([ s *' s i+iL K d ). 

The space Conf is the space of configurations for the snakes of length L 
which are continuous and "piecewise C 1 for V" ■ The snake S z associated to 
z S Conf is the map S z (s) = Jqz(t)cIt. Taking its "snout" S Z (L) provides a 
map / : Conf — > M. d , defined by 

f(z) = / z(s)ds 
Jo 

which is proven to be smooth. The image of / is the closed ball of radius L 
centered at the origin. 

The snakes corresponding to piecewise constant configurations, z(s) — z% for 
s G [sj_x, Si), are called polygonal snakes. In this case, we see / as a map from 
^gd-i-jjv to jjd serLC iiiig z = (zi, ■ ■ . , zjv) to J2^ =1 (si — Si-x)z{. For more details 
on this particular case see Ha2 . If all the (sj — Si_i) are equal, the snake is 
called an isosceles polygonal snake. 

The critical points of / are the lined configurations, where {z(s) \ s G 
[0,L]} C {±p} for a point p G S d_1 (for polygonal snakes, this is [Hall Theo- 
rem 3.1]). The snake associated to such a configuration is then contained in the 
line through p. The set of critical values is thus a finite collection of (d — 1)- 
spheres centered at the origin, which depends on V ■ 

Charming snakes will now be a path-lifting ability for the map /: given 
an initial configuration z G Conf and a C 1 -curve 7 : [0, 1] — > M d such that 
7(0) = f{z), we are looking for a curve ( h z, £ Conf such that zq = z and 
f(z t ) = 7(£). In Ehresmann's spirit, we are looking for a connection for the map 
/. The tangent space T z Conf to Conf at z is the vector space of those maps 
v G C^([0,L],R d ) such that (v(s),z(s)) = for all s G [0,L], where (,) denotes 
the usual scalar product in M. d . We endow T 2 Conf with the scalar product 
(v, w) — j^(v(s),w(s))ds. For each smooth map ip : R d — > K, one gets a vector 
field Grad(t/?o/) on Conf defined by 

Grad z ((po/)(s) =grad /(z) ^- (z(s), grad f(z) <p)z(s) . 

This vector field plays the role of the gradient of tpo f, that is 

(Grad^o /),!,) =T s (<pof)(v) 

for each v G T z Conf (as the metric induced by our scalar product is not com- 
plete, gradients do not exist in general). For z G Conf , the set of all Grad z (<p° f) 
for ip G C 1 (R d , R) is a vector subspace A z of T z Conf , of dimension d — 1 if z is 
a lined configuration and d otherwise. The correspondence z t— > A z is a distri- 
bution A (of non-constant dimension). For a pair (z,j) G Conf x C 1 ([0, 1],R ) 
such that f(z) — 7(0), the snake charmer algorithm takes for z t the horizontal 
lifting of 7 for the connection A. 
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As the map / is not proper and the dimension of A z is not constant, the 
existence of horizontal liftings has to be established. We use the C°°-action of 
the Mobius group M6b(d — 1) on Conf by post-composition: g ■ z = goz. For 
z £ Conf , let A(zq) be the subspace of those z £ Conf which can be joined 
to z by a succession of A- horizontal curves. One of the main results f jRol 
Theorem 2.19], proven in |Ha2| for isosceles polygonal snakes) says that A(zq) 
coincides with the orbit of zq under the action of M6b(<i — 1): 

Theorem 2.2 For all Zq £ Conf, one has A(zo) — M6b(d — 1) • Zq. 

The proof of Theorem 12.21 uses the following fact about the action of the 1- 
parameter subgroup 17 of M6b(e? — 1) introduced in II. 21 

Lemma 2.3 The flow <fit(z) — ■ z on Conf is the gradient flow of the map 

z ^ (/(«))«>• 

The flow <f>t on Conf is therefore A-horizontal. If z £ Conf , denote by f3 z : 
Mob(d — 1) — > Conf the smooth map /3 z (g) — g ■ z. Lemma [2.31 implies that 
the map f3 z sends the distribution A H onto the distribution A. If z is not a 
lined configuration, T g j3 z : A^ — > A g . z is an isomorphism; if z is lined, there is 
a 1-dimensional kernel (for instance, if z(s) £ {±p}, then Y\ ■ z = z for all t). 

2.4 The differential equation. As a consequence of Theorem l2.2l a A-horizontal 
curve Zt starting at z$ may be written as g(t) ■ zq where g(t) G M6b(c?— 1) and 
g(0) = id. For a given C 1 -curve j(t) £ M. d , the A-horizontal lifting g(t)-Zo of 7(4) 
is obtained by taking for g(t) the solution of an ordinary differential equation 
in the Mobius group M6b(d — 1) that we describe here below (more details can 
be found in |Rol Chapter 3]). 

For z £ Conf, we define the (d x d)-matrix M(z) by 

7 L z 1 (s)z 1 (s)ds ■■■ J o L z 1 (s)z d (s)ds\ 

,/ zi(s)z d (s)ds ■■■ J z d (s)z d (s)dsj 

(Observe that the second term is the Gram matrix of the vectors {z\, . . . , z d ) 
for the scalar product). It can be proved that M(z) is invertible if and only if 
z is not lined. Let F = F ZQ : M6b(d — 1) — > R d be the composition F = f°/3 Zo . 
Let g £ M6b(d— 1). Using the linear injective map \ of ll-21 with image Tt, and 
the right translation diffeomorphism R g in M6b(ii — 1), we get a linear map 

R d -4- H ^ A« ^ T F(g) R d « R d (2) 

It turns out that the matrix of the linear map 1S M(g ■ zq). If g ■ zq is not 
a lined configuration, T g F is bijective and M(g ■ zq) is invertible. One way to 
insure that g ■ zq is not lined for all g £ M6b(rf — 1) is to assume that zo takes 
at least 3 distinct values. The following result is proven in [Rol Prop. 3.10]: 
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Proposition 2.5 Let zq £ Conf be a configuration that takes at least 3 distinct 
values. Let 7 : [0, 1] -> R d be a C 1 -curve with 7(0) = f(z). Suppose that the 
C 1 -curve g : [0, 1] — > M6b(d — 1) satisfies the following differential equation: 

g(t)=R g(t) o X (M- 1 (g(t)-z )i(t)) , 5 (0)= id. (3) 

Then, the curve g(t) ■ zq is the unique /^.-horizontal lifting of 7 starting at Zq. 

As the Mobius group is not compact and the map F is not proper, Differ- 
ential equation © may not have a solution for all t € [0,1]. The notion of 
"sedentariness", described in 12. 71 is the main tool to study the global existence 
of A-horizontal liftings. 

2.6 The cases d — 2 or d = 3. For planar snakes (d = 2), one can use 
that Mob(l) is isomorphic to P5T7(1, 1) = SU(1, l)/{±id}, where 517(1,1) = 
{(5 a) I a ' & e C and H 2 ~ \ b \ 2 = The isomorphism between M6b(l) and 
PSU(1, 1) comes from the SU(1, l)-action on C by nomographic transforma- 
tions, which preserves the unit circle. The group 577(1,1) being thus a 2-fold 
covering of M6b(l), the curve git) £ M6b(l) of Proposition^}] admits a unique 
lifting in g(t) £ SU(1, 1) with g(0) — id. The output of the snake charmer 
algorithm will then be of the form z% = g(t) ■ z , using the action of SU(1, 1) 
on Conf (1). Working in the matrix group 577(1, 1) is of course favorable for 
numerical computations. 

The Lie algebra psu(l, 1) = su(l, 1) consists in matrices of the form ( ™ _ iu ) 
with a£l and b £ C. Under the isomorphism Slt(l, 1) « mdfa(l), the element 
x(v) £ H C mob(l) corresponds to the matrix (Q g) £ su(l, 1). These matrices 
form a 2-dimensional vector space H. C su(l, 1) giving rise to a right invariant 
distribution A n on SU(1, 1). By Proposition 12. 51 the curve g(t) ■ zq £ Conf (1) 
is horizontal if g(t) is A w -horizontal. In this language, Equation © becomes 

where u = v%(g,t) + 1^2(5, i) is obtained by 

/ Wl ( ff (t),i)\ _ M ,- m 2 ^7i(*)\ 

For more details, see |Eol Chapter 5]. An analogous (and formally similar) 
matrix approach is available for spatial snakes (d = 3), using the isomorphism 
PSL(2,C) m Mob(2); see (H5I Prop. 5.6]. 

2.7 Sedentariness. Denote by /1 the Lebesgue measure on M. Let z € Conf . The 
sedentariness sed(z) £ [0,L] of z is defined by scd(z) = max peSt i-i /i(z _1 (p)). 
This maximum exists since the set {p £ § d_1 | /i(z _1 (p)) > r} is finite for 
all r > 0. By Theorem 12.21 sed(z) is an invariant of A{z). Observe that, if 
sed(2:) ^ -7/2, then z takes at least 3 distinct values. The sedentariness of z is 
used in |Rol Section 3.3] to get global existence results for horizontal liftings of 
a path starting at f(z). The one we need is the following 
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Proposition 2.8 Let z £ Conf and let 7 : [0, 1] — > R d be a C x -path with 7(0) = 
f(z). Suppose that j([0, 1]) is contained in the open ball centered at of radius 
L— 2sed(z). Then there exists a (unique) ^-horizontal lifting 7 : [0, 1] — ► Conf 
for 7 with 7(0) = z. 

As the sedentariness of a configuration is preserved along horizontal curves, 
Proposition 12. 81 is also true for continuous piecewise C 1 -paths. A configuration 
z is called nomadic if sed(z) = 0. Proposition 12.81 guarantees that, if z is 
nomadic, any C 1 -path starting at f(z) admits a horizontal lifting, provided its 
image stays in the open ball of radius L. More general results may be found 
in |Rol Section 3.3]. 

2.9 Continuity of the algorithm. We now describe how the snake charmer algo- 
rithm behaves as we vary the initial configuration zq and the C '■-curve 7. Our 
goal is not to present the most general statements but only those needed in 
Section |21 

Let zq £ Conf, b = f(zo) and let Conf b = Set a = sed(zo) and 

suppose that |/(«o)| < L — 2a. Denote by _Bl_ 2(T (0) the open ball in R d of 
radius L — 2a centered at the origin. Consider a C 1 -curve 7 : [0, 1] — > Bl-2<t(0) 
with 7(0) = b. Define E a j, = {z 6 Conf b | sed(z) < a}; it is a metric space 
with the induced metric from Conf . Consider the map 

X : M6b(d - 1) x E a>b x [0, 1] 1 — > TM6b(d - 1), 

defined by 

X(g,z,t) - Rg^ox^Af-^git) ■ z ) 7O)) • 

This is a vector field on M6b(d— 1) depending on the time t and on the parameter 
z e E a £. Differential equation J3J becomes 

g(t) = X(g(t),z,t) , 5 (0) -id. (5) 

For any given z E E a j,, the solution g z of Equation J5J exists for all t £ [0, 1] by 
Proposition 12.81 Since X is continuous and its derivative in g is a continuous 
map on the variables z and t, classical results on the dependence of parameters 
for the solution (see for example |Scl 4.3.11]) imply that the map (z, t) ^ gz(t) 
is continuous. This yields: 

Proposition 2.10 (Continuity in z) Let 7 : [0,1] — > £?l_2o-(0) be aC x -curve 
such thatj(0) = b. For any z £ E a ^, denote byj z (t) = g z {t)-z the ^-horizontal 
lift 0/7 starting at z. The map E G j, x [0, 1] — * Conf that sends {z,t) to 7*(t) 
is continuous. 

Let M C E ai b be a smooth submanifold of Conf and suppose that 7 is of 
class C 2 . The map X restricted to M is therefore C 1 (in all variables) and using 
once again jSc] 4.3.11] we get: 
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Proposition 2.11 (Differentiability in z) Let 7 : [0, 1] — > -E>l-2o-(0) &e a 
C 2 -curve such that 7(0) = b. For any z £ M, denote by 7 2 (t) = fjz(t) ■ z the 
^.-horizontal lift 0/7 starting at z. The map M x [0, 1] — > Conf t/iat sends (z, t) 
to %(t) is of class C . 

We have an analog result when we fix the configuration z$ and vary 7. 
Consider the set C a , b = C^O, l],M d ) | 7 (0) = b, 7 ([0,1]) C B l _ 2ct (0)} 

that is an open subset of the affine space of C 1 -paths starting at b. 

Proposition 2.12 (Continuity in 7) Let zq £ Conf , b = f(zo) and a = 
sed(zo)- Suppose that |/(zo)| < L — 2a. For any 7 £ C a f,, denote by 7z (t) the 
/^.-horizontal lift of 7 starting at zq. The map Ccf, x [0,1] — > Conf that sends 
(7, t) to j z (t) is continuous. 

2.13 Bivalued configurations. Let z be a bivalued configuration, that is zo([0, L]) = 
{po, go}, with pq =/= qq. Let L p and I/ g be the Lebesgue measures of Zq 1 (po) and 
Zq (go) (these preimages are finite unions of intervals). By Theorem 12 .21 a hor- 
izontal curve Zt starting at zq will stay bivalued: z t {[Q,L\) = {p(t),q(t)}, with 
p(0) = po, g(0) = g and p(t) ^ q(t). Also, one has z t _1 (p(t)) = Zq 1 {p ), 
z t _1 (g(i)) = (9o) and therefore zt is determined by the pair (p(t),q(t)). 
Hence, A(zq) is contained in a compact submanifold W of Conf naturally 
parameterized by § d_1 x § d_1 . The restriction of / to W gives a C°°-map 
/ : W — > R d which takes the explicit form /(p, g) = L p p + i g g. The definition 
of the distribution A implies that A( p iQ \ C Tq^W when (p, q) £ W. 
Consider the open sets W° C W and U° C K d defined by 

W° :={(p,g) eW\p^ ±g} , C/° := {1 £ l d \x\ < L and \x\ ^ L p - L q } 

As W° contains no lined configurations, the map / restricts to a submersion 
/ = / : W a — > U for which A is an ordinary connection. As / extends 
to / : W — > R d and as W is compact, the map f° is proper. The original 
Ehresmann's construction of horizontal liftings |Eh| then apply: if z £ W° and 
7 : [0, 1] — * U° is a C 1 -curve 7(0) = zq(L), then 7 admits a A-horizontal lifting 
7 : [0, 1] -> of /° with 7(0) = z . 

Let us specialize to planar snakes (d = 2). Each fiber of /° then consists of 
two points; as f° is proper, it is thus a 2-fold cover of U°. Therefore, 7 (t) is 
determined by /(7(f)) = 7(0- We deduce that t/ie unique lifting 0/7 into Conf 
is horizontal. 

If, for f = to, 7(t) crosses the sphere of radius L p — L g , it may happen 
that 7(t) tends to a lined configuration when t — > to (see Example 14.21 below). 
To understand when the unique lifting 7(f) is horizontal at t$, we must study 
horizontal liftings around a lined configuration, which, after changing notations, 
we call again z$, corresponding to (j>q, go) £ W with p Q = — q . The vector space 
A Zo is then of dimension 1 and it turns out that T Zo f(A Zo ) is the line orthogonal 
topo, see |Hq1 Remark 1.17]. Let 7 : [0,1] -> R d be a C^curve with 7(0) = z (L) 
and 7(0) 7^ 0. A necessary condition for 7 to admit a horizontal lifting is then 
(7(0), po) = 0. If we orient the plane with the basis (po,7(0)), the curve 7 
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has a signed curvature k(0) at t = 0. In |Rol Prop. 3.18], it is proven that, 
around t — 0, the unique lifting of 7 into Conf is horizontal if and only if 7(0) 
is orthogonal to po and 

K(0 ) = hi_hiL. (6) 

Condition © has been detected by the numerician Ernst Hairer. Such a second- 
order condition for the existence of a horizontal lifting for a non-constant rank 
distribution is worth being studied. As far as we know, no such a phenomenon 
is mentioned in the literature. 

2.14 Miscellaneous. We finish this section by listing a few more properties of 
the snake charmer algorithm. Let (5, 7) be an input for the algorithm, with S a 
snake of length L. Let St ■ [0, L\ — > M d be the output. The following two results 
follow directly from Theorem 12. 21 

Proposition 2.15 (Regularity) Suppose that, on some open subset U C [0,L], 
the snake S is of class C k , k £ {1,2,..., 00}. Then the same holds true for the 
snake St for all t € [0, 1]. 

Proposition 2.16 (Periodicity) Suppose that there exists T £ M. such that 
z(s) = z(s + T) for all s such that s and s + T belong to [0, L]. Then, Zt{s) = 
z t {s + T) for all te [0,1]. 

The following proposition follows either from Theorem 12.21 or simply from 
the fact that St is a horizontal lifting. 

Proposition 2.17 (Reparameterization) Let ip : [0, 1] — * [0, 1] be an orien- 
tation preserving C 1 -diffeomorphism. Then, the deformation of S following the 
curve j(ip(u)) is S^ 

Finally, by construction, the distribution A is orthogonal to fibers of / : 
Conf -> R d . This can be rephrased in the following 

Proposition 2.18 z t is the unique lifting of 7 which, for all t, minimizes the 
infinitesimal kinetic energy of the hodograph, that is 5 J^\-^Zt(s)\ 2 ds. 

Again, the existence of such minimizer St is only guaranteed by an analysis 
like inl2"7l 



3 Holonomy orbits 

Let b £ M d . Define Conf b — the space of all configurations with as- 

sociated snake ending at b. Define the holonomy orbit horb(zo) C Conf b of 
zq e Conf b by horb(zo) = A(zq) n Conf b ; it is thus the subspace of those 
z £ Conf which are the result of the holonomy of the snake charmer algorithm 
for a pair (zq,j) with 7 a piecewise CMoop at b. Even if Zq is not lined, the 
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point b might not be a regular value of /, so might not be a submani- 

fold of Conf . But Lemma f3. II below tells us that this is the case for horb(zo). 
For z € Conf , define the spherical dimension spdim(z) of z to be the minimal 
dimension of a sub-sphere of containing the set z([0, L\). By Theorem 12 .21 
spdim(z) is an invariant of A(z). Notice that spdim(z) = if and only if the 
configuration z takes one or two values. 

Lemma 3.1 Let zq G Conf . Let k = spdim(zo) and suppose that k > 0. Then 
horb(zo) is a smooth submanifold of Conf of dimension ^] (d — i). 

Generically, spdim(zo) = d— 1, so horb(zo) is of dimension d(d — l)/2. The case 
spdim(zo) = is not covered by Lemma l3.ll It contains the monovalued case (z 
being constant) where horb(z ) = f~ 1 (f( z o)) — { z o}- The other case, formed 
by the bivalued configurations, is interesting and is treated in Example 14.21 and 
Proposition 14.41 

Proof of Lemma Em Let (3 : M6b(d— 1) — > Conf be the map /3(g) = g ■ z 
and let F : Mob(d-l) -> R d be the composition F = fo/3. Let K = F- 1 (f(z )). 
By Theorem 12 .21 one has horb(zo) = /3(K). The map F is smooth. 

Since spdim(zo) > 0, the configuration zq takes at least three values in S d_1 
and so does g ■ zq for all g g M6b(d — 1). Therefore, A(zq) contains no lined 
configurations and the tangent map T g F : T g M6b(d - 1) T F ( g) R d w R d is 
surjective for all g € M6b(d— 1) (restricted to A^, it is an isomorphism). Hence, 
F is a submersion. Thus, if is a smooth submanifold of M6b(<i— 1) of dimension 
dim M6b(d - 1) - d = d(d + l)/2 - d = d(d - 1) /2. 

The manifold K contains the stabilizer A of zq and K ■ A — K . Let V k be 
the smallest sub-sphere of S d_1 containing the image of Zq. The group A is then 
the stabilizer of the points of V k ; since k > 0, A is conjugate in M6b(rf — 1) to 
SO(d — k — 1). Therefore, the quotient space K/A is of dimension 

i=l j=l i=l 

By |Rol Proposition 2.33], the map (3 induces an embedding of M6b(c? — l)/A 
into Conf, hence an embedding of K/A into Conf with image horb(zo). This 
proves Lemma l3~Tl [] 

In general, horb(z) is not closed; see examples in |Ha2l pp. 1 12-1 13] 1 . But 
this is the case if f(z) is near the origin: 

Theorem 3.2 Let zq £ Conf such that \f(zo)\ < L — 2sed(zo)- Let k — 
spdim(zo) and suppose that k > 0. Then, each connected component o/horb(zo) 
is a closed smooth submanifold of Conf which is diffeomorphic to the homoge- 
neous space SO(d) / SO(d — k— 1). 

x In |Ha2) . one studies the set Hb(U) ■ zo: if 6 = /(zo) and U is an open neighborhood of 
b, 1i. b (U) denotes the subgroup of diffeomorphisms of Conf b obtained as the holonomy of a 
piecewise CMoop at b that stays in U. Notice that 7it(U) ■ zq C horb(zo) and that Hf,(U) ■ zo 
contains the connected component of zo of horb(zo). 
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Observe that SO{d) / SO(d — k — 1) is the Stiefel manifold of orthonormal 
(k + l)-frames in R d . 

Proof: Suppose first that /(zo) = 0. Let z G horb(z ). For each g e SO(d), 
one has f(g ■ z) — g ■ f(z) = 0, hence SO(d) ■ z C horb(zo) = horb(z). As in 
the proof of Lemma l3~Tl let [i z : M6b((i — 1) — > Conf be the map f3 z {g) = 9 • z. 
As spdim(z) = spdim(zo) = k > 0, the stabilizer of z is conjugate in SO(d) to 
SO(d— k — 1). Therefore, (3 z (SO{d)) is a connected, compact submanifold of 
Conf diffeomorphic to SO{d)/SO(d -k-1). By Lemma EH1 (3 z (SO(d)) must 
be the connected component of z in horb(zo)- This proves Theorem 13.21 when 
f(z ) = 0. 

In the general case, let b = f{zo) and let 7 : [0, 1] — » R d be the linear 
parameterization 7(f) = tb of the segment joining to b. li 5 : [0, 1] — > X 
is a map, one denotes by <5 - : [0,1] — > A" the map <5~(t) = 5(1 — t). Let 
2 G Conf b such that sed(z) < sed(zo)- By Proposition ^. 81 the path 7 admits a 
A-horizontal lifting j z in Conf, starting at z. The A-parallel transport t(z) — 
7 Z (1) is then defined. Of course, 7" also admits a horizontal lifting 7~( z ) which 
is equal to (j z )~ ■ By Proposition 12.101 the parallel transport gives rise to a 
homeomorphism 

r : {z e Conf 6 | sed(z) < sed(z )} Conf | sed(z) < sed(z )} . 

One has r(horb(zo)) = horb(r(zo)). As horb(zo) is a smooth submanifold of 
Conf, 

i>orb(z ) &n( ^ r |r(horb(zo)) are smootn by Proposition 12. 1 II Therefore r 
induces a diffeomorphism from horb(zo) onto horb(r(zo))- This terminates the 
proof of Theorem 13. 21 [] 

In the particular case of planar snakes (d — 2), we get: 

Corollary 3.3 Let zq S Conf (1) such that |/(zo)| < L — 2sed(zo) and 
spdim(z ) = 1. Then horb(zo) is a disjoint union of circles. 

We do not know, under the hypotheses of Theorem 13.21 or Corollary 13.31 
whether the manifold horb(zo) is connected. This is not true if |/(^o)| > L — 
2 sed(zo), see Example 14.21 It is however the case if sed(zo) = 0: 

Proposition 3.4 Let zq € Conf be a nomadic configuration. Letk = spdim(zo)- 
Then, horb(zo) is a closed smooth submanifold o/Conf which is diffeomorphic 
to the homogeneous space SO(d)/ SO(d — k — 1). 

To prove Proposition l3.4l we need the following lemma whose proof is post- 
poned till the end of this section. 

Lemma 3.5 Let z$ G Conf and let z £ horb(zo)- Then, there exists a loop 
7 : [0, 1] — > M. d at /(^o) which is of class C°° and which admits a horizontal 
lifting joining z$ to z. 
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Proof of PROPOSITION 13.41 A nomadic configuration is not lined. There- 
fore, k > and the condition |/(^o)| < L — 2sed(z ) = L is automatic. By 
Theorem 13.21 it is then enough to prove that horb(zo) is connected. 

Let z S horb(zo). By Lemma 1531 there exists a loop 7 : [0, 1] — > M. d at f(zo), 
of class C°°, admitting a horizontal lifting joining z$ to z. For s S [0, 1], define 
7a (t) = (1 — 5)7(0) + 57(4). The map s 1— > 7 S is a homotopy of C°°-loops at 
f(zo), from 7 to the constant loop. By Proposition 12.81 each loop 7 S admits a 
horizontal lifting -% starting at Zq. By Proposition 12. 12*1 the curve s <— > 7«(1) 
is continuous, producing a path in horb(zo) from z to Zq. This proves that 
horb(zo) is connected, provided we prove Lemma [331 [] 

PROOF of Lemma "331 By Theorem l2~2l there exists g £ M6b(rf - 1) with 
z = g ■ Zq. The Mobius group M6b(d — 1) is arc-wise connected using only 
curves that are piecewise trajectories of flows of the type J?j (see |Ha2l p. 102]). 
Therefore, there exists (v u h) € M. d X E, i = 1, . . . , r, such that g = • • • . 
For £ = 1, ... , r, let & : [0, 1] [0, 1] be a C°°-function such that ^([0, ^-]) = 
{0} and 1]) = {1}. The map g(t) : [0, 1] -> M6b(d - 1) given by 

= F ^(t)A r ' '• r ^ 1 (t)A 1 

is a C°°-curve in M6b(<i — 1), joining the unit element to g. By Lemma 12.31 
the curve g{t) is A w -horizontal. Then g(t) ■ zq is the horizontal lifting of the 
C°°-path -f(t) = f(g(t) ■ zq), which is a loop of class C°° at f(zo). fj 

4 Examples 

All our examples consist of planar snakes (d — 2). We identify M 2 with the 
complex plane C. Configurations z G Conf (1) are expressed under the form of 
s 1 > e te ^ with 6(s) £ M. 

4.1 In our first example, the snake is a half circle with configuration z : [0, n] — > 
S 1 given by z(s) = ie~ is , thus 5 = 5 Z is given by S(s) = 1 — er %s . We have 
f(z) — S(tt) — 2. The curve 7 is a small circle, centered at (2.1875,0) and 
of radius 0.1875, followed in the trigonometric direction. The snake charmer 
algorithm has been solved using the method of 12.61 After one turn, the snake 
slightly leans to the left. Figure A below shows the snake after various numbers 
of turns of 7. One sees that after 326 turns, the snake seems being back in its 
initial position. 



11 




Figure A 

The curve t i— ► g(t) <E M6b(l) obtained from Equation © (such that z t = 
g(t) ■ zq), may be visualized using the diffeomorphism M 2 x S 1 M6b(l) 
described at the end ofO The two fold covering SU (1,1) -> K 2 x S 1 thus 
obtained is given in formula by 

with = 2arg(a) and v = 2 arccosh(|a|)e iars '- ab - ) . Figure B below illustrates the 
subset {(v(n ■ 2ir),9(n ■ 2n)) \ n= 1, . . . , 326} C K 2 x S 1 . The picture on the 
left hand side shows the 326 points of the set {v(n ■ 2ir)} (v = (ui,^)) and 
the one on the right hand side is the graph of the function n i— ► 9(n ■ 2tt) (also 
formed by 326 points, but hardly distinguishable). That the points look more 
concentrated between 80 and 250 seems to be related to the fact, seen in Figure 
A, that the snake's shape changes less drasticly in this range. As in Figure A, 
we observe that 9(t) and v(t) seem to have returned to their original position 
after 326 turns. 
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V2 



t v(0) 



3 




Figure C below shows the entire curve (6{t),v(t)) for t G [0, 326 • 2ir}. 




Figure C 

4.2 A non-connected holonomy orbit. Let z : [0, 2y/2] — > S 1 be the 2-valued 
configuration 

e i7r/4 if s < V2 

e-W 4 if S > . 



zo(s) 



Thus, /(z) = 2. Let 7 : [0, 1] -> C be the piecewise C°°-loop at 2 defined by 



7(*) 



2e 2,rt if t < 1/2 
8(i-3/4) if t > 1/2 
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By 12.131 the unique possible lifting z t , starting at zq and determined by 
f(z t ) = 7(f) is horizontal. It follows that zi = zq (the complex conjugate of zq). 
Hence horb(zo) = {zo,zo} is not connected. 

Using an arbitrarily small bump around 0, one can perturb 7 so that it avoids 
the origin. But z\ would then be equal to zq. So z t does not depend continuously 
on 7, contrarily to what happens in Proposition 12.1 21 Here, two hypotheses of 
Proposition ^. 121 are not satisfied: spdim(zo) = and, what seems more serious, 
the condition |/(zo)| < L — 2sed(zo) is impossible since sed(zo) = y/2 — L/2. 

4.3 Example l4.2l is an illustration of the following computation of horb(zo) when 
zo is bivalued. 

Proposition 4.4 Let zq be bivalued configuration and let b — f(zo). Then 
(i) horb(zo) ~ E> d ~ 2 if Zq is not lined; 
(ii) horb(zo) is one point if Zq is lined and f(zo) =/= 0; 
(Hi) horb(zo) ~ S^ 1 if Zq is lined and f(zo) = 0. 

Proof: Let p , q , L p , L q , f : W -> R d and W° as in |2~T31 One has 
horb(zo) C /^(fo). If w E R d , denote by B w be the stabilizer of w in SO(d). 
We divide the proof into four cases. 

Case (i) with d = 2. Let 7 be a loop at b, with |7(f)| < L, such that 7 hits 
the sphere of radius \L p — L q \ in a single point, tangentially and respecting the 
curvature formula © . Bv l2.13l 7 has a A-horizontal lifting z t in W and Zt inter- 
sects the submanifold W° of lined configurations in one point and transversally. 
This implies that z\ ^ zq, which proves Case (i) when d = 2. Example 14.21 
illustrates this argument with L p = L q . 

Case (i) with d > 2. The vectors p, q and b are co-planar and p ^ b ^ 0, so 
Bf, ■ zq f=a B b /(B b n B p ) w SO(d - l)/SO(d - 2) w § d - 2 . One has B b ■ z C 
horb(zo) C / _1 (6). Case (i) with d > 2 then follows from the fact that / _1 (6) w 
gd-2 j a ^ er j g p r0 ven in jl la 1 Prop. 4.1] but, for the convenience of the 

reader, we repeat the argument here in our language. Let W b = / _1 (lR>o&), 
which is a submanifold of W of dimension d— 1 (see Hal , (1.3)]). Let h : W b — > R 
defined by h(z) = |/(z)|, which is a Morse function whose critical point are the 
lined configurations of W b (see |Hall Theorem 3.2]). If z e W b is a critical point 



14 



of h with h(z) > b, then z — (b/\b\,b/\b\), a non-degenerate maximum, and 
h(z) = L > |6|. Therefore, = h^Qb]) is diffeomorphic to § d ~ 2 by the 

Morse Lemma. 

Case (ii). This case is trivial since consists of the single point (po, —po)- 

Case (Hi). The map § rf_1 — ► W given by p i— ► (p, — p) is an embedding with 
image / _1 (0), as well as the inclusion B ■ z c horb(z ) C j _1 (0). [] 
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